Astronomy & Astrophysics manuscript no. 19422 


© ESO 2013 


March 6, 2013 





The size, luminosity and motion of the extreme carbon star 

IRC+10216 (CW Leonis) 

K. M. Menten 1 , M. J. Reid 2 , T. Kamiriski 1 , and M. J. Claussen 3 



(N 

o 

S3 
33 

o 

m 

6 



(N 
> 
OS 

m 

in 
o 

(N 



X 



1 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany e-mail: kmenten, 
kaminski@mpif r-bonn . mpg . de 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street/MS42, Cambridge MA 02138, USA e-mail: 
reidOcf a . harvard . edu 

3 National Radio Astronomy Observatory, Array Operations Center, P.O. Box O, Socorro, NM 87801, USA e-mail: 
mclausse@nrao . edu 



Received / Accepted 



ABSTRACT 



Very Large Array observations of the extreme carbon star IRC+10216 at 7 mm wavelength with 40 milli-arcsecond 
resolution resolve the object's radio emission, which forms an almost round uniform disk of 83 milli arcseconds diameter, 
corresponding to 11 AU (for an assumed distance of 130 pc). We find a brightness temperature of 1630 K for the radio 
photosphere. Since the emission is optically thick, we can directly estimate IRC+10216's average luminosity, which is 
8600 1/0. This value is in excellent agreement with what is predicted from the period- luminosity relation for carbon- 
rich Miras. Assuming an effective temperature of 2750 K for IRC+10216, it implies an optical photospheric diameter 
of 3.8 AU. Our precise determination of IRC+10216's proper motion fits the picture presented by far-ultraviolet and 
far- infrared wavelength observations of its interaction region with the interstellar medium (its "astrosphere" ) : the star 
moves roughly in the direction expected from the morphology of the termination shock and its astrotail structures. 
Calculation of its three dimensional velocity and an analysis of the kinematics of its surrounding interstellar medium 
(ISM) suggest an appreciable relative velocity of 42 km s _1 , which is about half the value discussed in recent studies. 
This suggests a lower (time-averaged) mass loss rate and/or a higher ISM density than previously assumed. 

Key words. ISM: molecules - Stars: circumstellar matter 



1. Introduction 

The carbon-rich evolved star IRC+10216 (also known as 
CW Leonis) is one of the most prominent an d best-studied 
near- infrared (NIR) sources in the sky (Bccklin et al. 
1969). IRC+10216 is probably a typical carbon star 
(|Herbie fe Zappalal [T97a iMillerl 11970ft ' near the end of its 
lifetime on the asymptotic giant bra nch (AGB), which 
is ch aracterized by extreme mass-loss (jHabing fc Olofssonl 



This and its proximity, make IRC+10216 a unique 
object of interest, allowing studies that would be very 
difficult or impossible for any other source of its kind. 
Estimated distances, D, inferred from modeling the 
CO emission in its envelope, range from 11 0-150 pc 
(|Crosas fc Mentenlll997HGroenewegen et al.lll998l ): we shall 
adopt a median value of 130 pc. 

IRC+10216's h i gh m ass-loss rate of 2 • 10~ 5 M yr" 1 
(jCrosas fc Menten] Il997l scaled to 130 pc) results in a 
dense circumstellar envelope (CSE) whose exceedingly 
rich chemistry can be easily studied at infrared 
through submill i meter t o rad i o wavelengths (se e, e.g. 
iCernicharo et~aH fl996l l2000t iPatel et alJ l20ll . The 
innermost part of the envelope, i.e., within ~50 AU or 20 
stellar radii, and the star itself is hardly observable at 
visual wavelengths, due to the high extinction of the dust 
that is abundantly produced in this region. In the infrared 
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(IR) regime, however, IRC+10216 presents a complex 
and dynamical picture. High spatial resolution imaging 
obtained with speckle intcrfcrometry and/or adaptive 
optics shows several distinct features on subarcsecond scales 
that vary over time-s cales of years, not only in position 
but a l so in luminosity (|Osterbart et alJl20fjd : IWeieelt et al.l 
2002; iLeao et aD l2006h . Combined data taken in the 
near/mid-IR H, K, L, M, and N bands, together with 
far-IR data, have been used to find a self-consistent m odel 
of the star and its envelope (Mcn'shchiko v et al.ll200ll ). 

All modeling efforts are hampered, however, by the 
inability to tell which (or if any) of the observed compact 
features actually corresponds to the stellar photosphere, 
with different approaches 

yiel ding widely differin g conclusions (see above references 
and lTuthill et al] |2005). M oreover, adaptiv e optics H-band 
imaging polarimetry by iMurakawa et al.1 (|2005l ) implies 
that the position of the illumination source (the central 
star) is different from any of the previously postulated 
positions. Needless to say, all radiative transfer modeling 
efforts of this keyst one envelope, sophist i cated as they 
may be (see, e.g. Ilvezic fc Elitzurl Il996t iGroenewegerj 
119971: iMen'shchikov et al.l l200lh severely suffer from the 
uncertainty of not knowing the position from which the 
luminosity originates. 

Clearly, a direct detection of the star that could 
unambiguously determine the stellar position would be 
superior to any of the above indirect methods. One of the 
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objectives of the present study is to determine IRC+10216's 
position at short radio wavelengths (7 mm) with an 
accuracy of a few milli arcseconds (mas) using the NRACQ 
Very Large Array. In a related project we seek to achieve 
absolute infrared astrometry to match our accurate radio 
astrometry. Since the radio position marks the bona fide 
location of the star, it will thus eventually be possible to 
unambiguously determine its IR counterpart. 

Radio emission from IR C+10216 has been studied at 
wave le ngths of 2 and 1 .5 cm (jSahai et al.lll989t iDrake et al.l 
119911) . iMenten et all (|2006t) presented Very Large Array 
observations at 3.6, 2, and 1.3 cm, which revealed an 
unresolved source < 95 mas and established that the 
emission is optically thick, i.e., its spectral index, a, is 
« 2; where flux density, S, is oc v a . This allows us to 
address a second goal of the present higher resolution 
(~ 40 mas) study, which is to determine the size and shape 
of IRC+10216's radio emission distribution and, using its 
implied brightness temperature, its luminosity. 

Our highly accurate position, combined with published 
data allow a high quality determination of IRC+10216's 
proper motion on the sky. This is an important quantity, 
given recent ultraviolet and far-infrared (FIR) observations 
of the interaction region between the star's expanding 
circumstellar envelope and the ambient interstellar 
medium. 

This paper has the following structure: In <J5] we give 
an account of our VLA observations. The determination 
of IRC+10216's size, accurate position, and proper motion 
are described in Sj3] The nature of its radio emission, 
the luminosity it implies, and its motion through the 
interstellar medium are discussed in SJU 

2. Observations and data processing 

2.1. VLA observations 

Our VLA observations took place on 2006 February 26 
(JD 2453793) with 23 antennas in operation. We used 
the maximum bandwidth setting provided by the VLA 
correlator. This comprises 2 intermediate frequency (IF) 
bands, each with 43 MHz effective bandwidth. We recorded, 
both, right and left circular polarization in each IF. One IF 
was centered at 43.3149 GHz, the other at a 50 MHz higher 
frequency. 

In order to ensure optimal calibration for our high 
frequency data, we employed the " fast switching" technique 
first described bv lLim et all (119981) . Over 7.5 h, we switched 
between IRC+10216 and one of the two nearby calibrators, 
J0943+170 or J0954+177. We used a cycle with 60 seconds 
on IRC+10216 and 40 seconds on one of the calibrators, 
which we repeated for 1-1.5 h segments. In order to check 
the effectiveness of our observing mode, between segments 
we switched between the two calibrator sources for ~ 10 
min. 

We improved the VLA antennas' pointing accuracy 
by updating pointing offsets once an hour, by observing 
J0954+177 in "interferometer pointing mode". For higher 
sensitivity these pointing scans were taken in X-band (8.4 
GHz). The corrections to the pointing model were applied 
real time. 

1 The National Radio Astronomy Observatory (NRAO) is 
operated by Associated Universities, Inc., under a cooperative 
agreement with the National Science Foundation. 



While processing the data we found that the position of 
J0943+170 as listed in the VLA calibrator database was 
insufficiently accurate, given the precision warranted by 
our data. We adopted better positions for both calibrator 
sources from the VLBI Global Solution 2008b Astro 
Catalog ma intained at the NA SA Goddard Space Flight 
CenteiQ; see lPetrov etaH (|2008l ). We used the task CLCOR 
within NRAO's Astronomical Image Processing System 
(AIPS) to correct our data for the improved positions. 
The adopted calibrator positions, as well as the position 
of IRC+10216 determined as described in are listed in 
Tabled) 

2.2. Calibration and initial uv-plane and image analysis 

Calibration and image processing was performed with 
AIPS. While loading the data into AIPS (using the 
task FILLM), corrections for the elevation-dependent gain 
curve and atmospheric transmission were applied. Visual 
inspection of the measured visibilities revealed a minimal 
amount of flawed data. The absolute flux density, S, scale 
was established by an observation of 3C286, which has 
S = 1.455 Jy at 43.317 GHz. Since this source is slightly 
resolved at our observing frequency, we used a clean- 
component model downloaded from the VLA website to 
determine amplitude and phase solutions using the AIPS 
task CALIB. Then we used CALIB to determine amplitude 
and phase solutions for J0954+177 and J0943+170, for 
which we bootstrapped flux densities of 0.234 and 0.106 
Jy, respectively, by comparing with the solutions found for 
3C286. Using CLCAL, we applied the amplitude and phase 
corrections to the IRC+10216 data. 

Using IMAGR we produced an image of IRC+10216 
and found the star offset in right ascension and declination 
direction, (Ax, Ay), by (+0."641, +0."157) relative to our 
phase-center position. To enable subsequent uw-domain 
analysis, we used UVFIX to shift the visibility data to the 
center of the icy-plane. We then used IMAGR to produce 
images of the shifted wv-dataset, both with uniform and 
natural weighting (see £)3.1[) 

3. Results 

3.1. Imaging of IRC+10216 

3.1.1. The size and brightness temperature of the radio 
emission 

Table shows the result of fits to IRC+10216's observed 
brightness distribution at 43.3 GHz. Fits were made to 
the image (using JMFIT) as well as to its corresponding 
ww-plane data (using OMFIT). We solved for the peak 
brightness and source size using elliptical models. Our 
results, represented in Tableland Fig. [21 indicate an almost 
circular source. Given this, we also performed fits with 
our independently developed least-squares fitting program, 
which assumes circular symmetry. This has the advantage 
that it delivers useful statistical quantities. An inspection 
of Fig. [5] indicates that the uniform disk model (blue line) 
matches the data much better than the Gaussian (red line). 
This is verified by a comparison of the reduced x 2 values, 
which are 1.35 for the disk and 3.07 for the Gaussian model. 
The best fit uniform disk has a diameter, d, of 83 ± 1 

2 http:/ /lacerta.gsfc. nasa.gov/vlbi/solutions/ 



2 



Menten et al.: Size, luminosity and motion of IRC+10216 
Table 1. IRC+10216 and calibrator sources 



Source 



a(J2000) 



<5(J2000) 



d.r 

n 



Ay 

(°) 



IRC+10216 (CW Leo) 

J0943+170 

J0954+177 

3C286 (J1331+3030) 



09 h 47 m 57. s 4255(6) 
09 43 17.22396(2) 
09 54 56.82362(5) 
13 31 08.28806(2) 



+13°16'43."815(10) 
+17 02 18.9628(12) 
+ 17 43 31.2222(2) 
+30 30 32.9592(5) 



-lfl 
+1?7 



S 

(Jy) 



0.013 
+3?8 0.106 
+4?4 0.234 

1.455 



Notes. Measured position for IRC+10216 (epoch 2006.16) and the adopted positions for the calibrator sources are given in columns 
two and three in J2000 coordinates. Calibrator positions are taken from the VLBI Gl obal Solution 2008b Astro Catalogue available 
on the NASA Goddard Spaceflight Center Geodetic VLBI group's website; see also lPetrov et all (|2008f ). Numbers in parentheses 
give the error in the last quoted digit(s). For the complex-gain calibrators J0943+170 and J0954+177, column four lists their 
angular separations from IRC+10216. The separation between these two calibrators is 2?9. Column six lists the measured flux 
densities based on the assumed value for 3C286. We estimate the flux density scale is accurate to ±5%. 



Table 2. Model fits to IRC+10216's 43.3 GHz brightness distribution 



Restoring beam Deconvolved source 



M #B,maj #B,min 


P.A.(B) S ,maj 


#S,min 


P.A.(S) 




/ SdO 


(mas) (mas) 


(°, E of N) (mas) 


(mas) 


(°, E of N) 


(mjy beam 1 ) 


(mJy) 


Imaging with natural weightin 


g - Gaussian model 










J 61 49 


-6.3 60 + 2 


54 + 2 


+34 ± 17 


6.22 + 0.07 


13.1 + 0.2 


13.1 + 0.2 2630 












Imaging with uniform weighting - Gaussian model 










J 44 39 


-6.8 63 + 2 


57 + 2 


+26±£ 


4.30 ± 0.09 


13.2 + 0.4 


13.3 ± 0.3 2450 










Visibility versus tti;-distance - 


Gaussian modeling 










O 


61 + 1 


55+3 


+31 ± 13 




12.8 + 0.2 


L 


56 ±2 








12.6 + 0.3 


Visibility versus itw-distance - 


uniform disk modeling 










O 


87 + 2 


'80 + 1 


+22 + 5 




12.2 + 0.1 


L 


83 + 1 








12.1 + 0.2 



Notes. The first four lines give results of Gaussian fitting, line 1 and 2 to a naturally weighted image and a uniformly weighted 
image, line 3 and 4 to visibility domain data. Lines 5 and 6 give the results for fits with a uniform disk model. The first column 
denotes the method (M) used: J for JMFIT and O for OMFIT (AIPS tasks). L stand for our independent least-squares program. 
For the first two lines, the second to forth columns give the FWHM major and minor axes and the position angle of the restoring 
beam determined by IMAGR, while the fifth to seventh columns list the same quantities for IRC+10216's brightness distribution 
and the eighth column gives fitted peak brightness. For all fits, the ninth gives the integrated flux density (the "zero spacing" flux). 
JMFIT and OMFIT fitted elliptical brightness models, whereas our least squares program used a circular model. 



mas, corresponding to 10.8 ±0.1 AU (D = 130 pc) and a 
brightness temperature, Tb, of 1635+82 K. The uncertainty 
in Tb is dominated by the 5% estimated uncertainty of our 
absolute calibration. 



3.2. Position and proper motion of IRC+10216 
3.2.1. Position 

Using JMFIT, we determined the best-fit position for 
IRC+10216 listed in Table[TJ The formal precision returned 
by JMFIT in the eastward, x, and northward, y, directions 
are 0.4 mas, much smaller than the values quoted in the 
table. To arrive at the listed, more realistic values for 
these quantities, we performed several tests. First, we only 
allowed phase and amplitude corrections determined for 
J0943+170 (in the following the "weak" calibrator) to 
be used to calibrate the IRC+10216 data. About half of 
the data, those taken while switching between IRC+10216 
and the "strong" calibrator, J0954+177, were discarded. 
Second, vice verse, we calibrated IRC+10216 only with the 
strong calibrator, discarding the data taken while switching 
with the weak calibrator. We then imaged both calibrated 



datasets and again used JMFIT to determine the resultant 
positions. We find position centroids with offsets, (Ax, Ay), 
relative to the position in Table[T]of (—2.8 mas, —6.4 mas) 
and (+2.8 mas, +5.9 mas), respectively. IRC+10216's best 
fit flux density determined from both of these sub-datasets 
agreed with that determined from the whole dataset to 
within 2%. 

To get a further estimate of the accuracy of our 
position determination, we calibrated the uv-data of the 
weak calibrator with the phase and amplitude corrections 
determined for the strong calibrator. Imaging of the former 
retrieved J09433+170's position with an offset of (Ax, Ay) 
= (—1.3 mas, —3.4 mas) relative to its nominal value 
quoted above. Calibrating, vice versa, strong with weak 
calibrator retrieved J0954+177's position at an offset of 
(+6.5 mas, +9.1 mas). In this exercise, both calibrator 
sources appeared unresolved and JMFIT returned formal 
upper limits of 7 and 10 mas for sizes of the weak and the 
strong calibrator, respectively. The flux densities of both 
were found to within a few percent of the bootstrapped 
value quoted above. We note that the weak and strong 
calibrator databases only contained 20% and 14% of the 
number of visibilities of the IRC+10216 database. 
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East Offset (arcsec) 



Fig. 1. Image of IRC+10216 at 43.3 GHz (6.9 mm). The 
brightness distribution is shown in grey scale. Contours 
give -5, 5, 10, 30, 50, 70, 90, and 99 percent of the 
peak brightness, which is 4.1 mJy beam -1 . The rms 
noise level in the image, 92 /zJy, corresponds to 2.2% of 
that peak value. The 41 mas FWHM symmetric restoring 
beam is represented in the lower left corner. The 83 mas 
diameter dashed circle gives the size of IRC+10216's radio 
photosphere derived by fitting a uniform disk brightness 
distribution to our VLA data. The dotted circle of diameter 
29 mas represents the diameter of the star's photosphere 
derived from the luminosity we determine and assuming 
an efficient effective temperature of 2850 K. Note that the 
major axes of Mars' and Jupiter's orbits (3.0 and 10.4 AU, 
respectively) would have angular diameters of 23 and 80 
mas at a distance of 130 pc. 

As listed in Table [TJ the separation on the sky between 
IRC+10216 and J0943+170 and J0954+177 is 3?9 and 5?9, 
respectively, and thus are 1.4 and 3.1 times larger than the 
2?9 arc between the two calibrators. Thus, it seems prudent 
to scale the position errors resulting from the calibrator- 
calibrator experiment by the ratio of the x and y differences. 
This would result in an absolute position uncertainty of ~ 3 
mas in each coordinate. We conservatively increase this to 
10 mas, because of the fact that our calibrators are both 
offset in declination and astrometric accuracy is expected to 
be worse for target-calibrator separations that are mostly 
in the north-south direction. 

3.2.2. Proper motion 

iBecklin et al.l (jl969fl . in the very first publication on 
IRC+10216, placed an upper limit of 30 mas yr™ 1 on 
its proper motion. These authors compared a position 
measured on a plate from 1969 with the position from 
an E-plate o f the Palomar Sky Survey taken in 1954. 
iMenten et al.l (|2006T ) used positions determined from VLA 
observations in 1987 and 1993 (with uncertainties of ss 25 
mas) to determine IRC+10216's proper motion in x- and y- 
direction, ([i x , /i y ), of (+26+6, +4+6) mas yr -1 . Combining 
these data with our more accurate position from 2006 
(Table 1), we calculate a proper motion of (+35+1, +12+1) 
mas yr -1 (see Fig. [3]). For a distance of 130 pc, these 




2000 4000 

uv — distance (kA) 



Fig. 2. Fringe visibility vs. baseline length for IRC+10216 
at 43.2 GHz (6.9 mm) on 2006 February 26. The red line 
represents the Fourier transform of a Gaussian brightness 
distribution with an FWHM of 56 mas and the blue line 
that of a uniform disk of diameter 83 mas (see text and 
Table©. 

angular motions correspond to linear speeds of 21.6 km s -1 
eastward and 7.4 km s -1 northward. This proper motion is 
a heliocentric value and needs to be corrected for the solar 
motion. The resulting values are 36.3 km s -1 eastward and 
13.2 km s -1 northward, implying a speed in the plane of 
the sky of 38.6 km s -1 along a position angle of 70° (E 
of N). We note that the motion we determine is in the 
opposite direction of that implied by t he binary model for 
IRC+10216, which [Guelin et all (|1993|) propose to explain 
an offset they find between the star and the centroid of 
the C4H emission imaged with the IRAM Plateau de Burc 

int erferometer . 

iRoeser et al.l (|2010l) combine data from the optical 
USNO-B1.0 all sky survey and the Two Micron All 
Sky Survey (2MASS) to calculate astrometric parameters 
for common sources, IRC+10216 amongst them. In the 
resulting PPMXL catalog, IRC+10216's proper motion is 
listed as (fi x ,fi y ) = (+28 ± 6,— 8 ± 6) mas yr -1 . Their 
eastward motion is in agreement with ours, while their 
northward motion is discrepant by w 3a. Given the 
complex structure of IRC+10216's innermost CSE and the 
significant variation of its morphology with wavelength, 
combined with the uncertainty of whether or not the stellar 
photosphere is observed at optical and IR wavelengths (see 
Sect. [1]), one might question the reliability of the PPMXL 
proper motion determination. 

4. Discussion 

4.1. Properties and nature of the radio emission 
4.1.1. Size and luminosity 

Our measured radio brightness temperature of 1660 K 
is definitely lower than IRC+10216's optical effective 
temperature, T». After critically evaluating models of 
AGB stars and a wide range of data (their Table 4) 
iMen' shchikov "et~aT1 (pOOl) conclude that T» lies between 
between 2500 and 2800 K. Thu s, our radio siz e and 
temperature are in line with what Rci d fc Menten! (|l997t 
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JD 2438000+ 

4000 8000 1.2 10 + 1.6 1 4 




1970 1980 1990 2000 2010 

Yea r 



Fig. 3. P roper motion of IRC +10216 determined for epochs 
1969.36 (iBecklin et al.lllQfiflh . 1987.51 and 1993.07 (both 
iMenten et al.ll2006t) . and 2006.16 (this paper). The red and 
blue dots mark the position offset relative to the epoch 
1987.51 position in the eastward (x) and northward (y) 
directions, respectively. The error bars for the 2006 position 
are comparable to the symbol size. Lines represents least- 
squares fits to the data. For the first two epochs the times 
have been slightly offset for clarity. 



their errors consist ent with the radio data discussed by 
Men ten et al.l (|2006l h Given this, a luminosity derived from 
radio data at any time is close to the star's average 
luminosity. Our estimate gains in reliability since we were 
able to determine IRC+10216's long-wavelength (infrared) 
phase, 0ir0, at the time of our VLA observations. As 
described in Appendix [21 we acquired data to determine a 
reliable value for fan, of 0.75. This means the luminosity we 
have determined is very close to the average value indeed. 

The (average) luminosity we determine, 8640 Lq (for 
D = 130 pc), is lower th an the value of 14000 Lq that 
iGroenewegen et ahl (|1998l ) find consistent with their multi- 
transition m odeling of the CO e missio n fro m IRC+10216's 
CSE while ICrosas k, Mentenl (|1997D and iDe Beck et all 
(|2012| ) find 11000 and 8500 L , respectively, in similar 
analyses (all for D = 130 pc). 

If we assume T cS = 2750 K for IRC+10216's 
effective temperature, i.e. , the median value considered by 
iMen'shchikov etlrtl (|2001h for their modeling, we calculate 
a diameter of 3.8 AU for the star's photosphere. This 
corresponds to f=s 1.3 times the major axis of Mars' orbit. 

The derived luminosity is comparable to the value one 
derives from the revised period- 

luminosity relation for carbon-ri ch Miras established by 
IGroenewegen fc Whitelockl (jl996lh For our best fit period, 
630 d (see Appendix |A"1). this relation predict s a luminosity 
of 983 Lq. Finally, we note that, remarkably. iBecklin et al.l 
(196§) in the very first publication on IRC+10216 derive 
10 31 W for its luminosity (assuming D — 200 pc), which 
corresponds to 11000 Lq for D = 130 pc. 



|2007|) find for the "radio photosphere," which surrounds 
the optical photospheres of M-type (oxygen-rich) AGB 
stars. They find that the radius of the radio photosphere 
is about twice the stellar radius and the temperatures is 
about l/v2 lower than at the stellar surface. Modeling 
the ionization balance and opacity sources in these r egions 
(mainly H free- free interactions). iReid fc Mentenl (|1997l ) 
predict a spectral index, anp, for the radio emission of 1.86. 

From measurem ents taken between 8 and 680 GHz, 
IMenten e t al. ( 200fJ) determined the radio-to-submillimeter 
spectral index of IRC+10216's radio emission to be 1.96 ± 
0.04, which is ste eper than the v al ue pr edicted for radio 
photospheres by IReid fc Mentenl (|1997|) and consistent 
with optically-thick blackbody emission in the Rayleigh- 
Jcans regime (cubb = 2). Since we are observing emission 
consistent with a blackbody, our measured diameter, d and 
brightness temperature, Tb, directly deliver IRC+10216's 
bolometric luminosity, L», via the Stefan-Boltzmann law: 
L* = nd 2 aT^, where a is the Stefan-Boltzmann constant. 
With T B = 1635 ± 82 K and d= 10.8 ± 0.1 AU (83 mas 
at D = 130 pc), we calculate L* — 3.3 • 10 37 erg s _1 , or 
8640 ± 430 Lq. The last uncertainty does not include a 
contribution from the distance uncertainty. 

Although IRC+10216's bolometric luminosity varies 
by a factor of 2.5 over a light cycle, model calculations 
predict that its radio flux varies by n o more than ±10% 
betw een minimum and maximum (|Men'shchikov et al.l 
l200l and pers. comm.). This is in excellent agreement 
with IRC+10216's 850 /jm light curve derived from data 
taken with the James-Clerk-Maxwell- Telescop43 and within 

3 http:/ /www .jach.hawaii.edu/JCMT/continuum/ 



4.1.2. Shape of the radio photosphere 

Performing VLA observations w ith angular res o lution 
and wavelength similar to ours, IReid fc Mentenl d2007h 
determined the shapes and sizes of the radio photospheres 
of three M-type AGB stars (o Ceti, R Leo and W Hya); 
the latter were all close to 6 AU. While o Ceti's radio 
photosphere has, like IRC+10216's, an almost circular 
shape, those of R Leo and W Hya appear to be significantly 
elongated. Formally, we find a flattening, e (= (a — b)/a) 
of 0.08 ± 0.02 for IRC+10216 (where a and b are the 
major and minor axis sizes, respectively). From Table 2 of 
IReid fc Mentenl (|1997ft we calculate e = 0.07 ± 0.13, 0.36 ± 
0.14 and 0.33 ± 0.14 for o Ceti, R Leo and W Hya, 
respectively 

IRC+10216's radio photosphere is almost perfectly 
round, comparable to that of the much lower mass-loss 
star o Ceti. Given its high mass-loss rate, IRC+10216 
has been conjectured to be on the verge of developing 
into a protoplanetary nebula (PPN) and asymmetries in 
its circumstellar material ha ve been interpreted a s signs 
for the onset of bipolarity (|Osterbart et al.l |2000T) . Here 
we emphasize that the radio emission is consistent with 
an extended photospheric origin and does not show the 
characteristics of a developing ultracompact HII region as 
observed in the more evolved PPN CRL 618. For this 



4 Conventionally, the phase, 4>, quoted for an AGB star is the 
time, measured as a fraction of its period, that has passed since 
optical maximum. Due to molecule formation, for most AGB 
stars, the IR maximum lags ~ 0.2 periods behind the optical 
maximum. This is well-esta blished for oxygen -rich Miras, but 
less so for carbon-rich ones ([Smith et al.ll2006l h 
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object, whic h has a ~ 30 times higher mass loss rate than 
IRC+10216 (|Young et alJl992ft . VLA observations show an 
enhanced radio flux associated with an extended ionized 
regio n that is elongated along the PPN's bipolar ou tflow 
axis (|Kwok fc Bignelllll984l: iMartin-Pintado et al.lll99l 

4.2. Motion through the interstellar medium 

Recent observations at ultraviolet, FIR and submillimeter 
wavelengths provide evidence for an interaction of 
IRC+10216's outflowing envelope and the ambient 
interstellar medium (ISM). Knowledge of the star's motion 
is an important parameter in any modeling of this 
ph enomenon. 

ISahai fc Chronopoulosl (|2010h present near and far 
ultraviolet (NUV and FUV) imagery obtained with the 
Galaxy Evolution Explorer (GALEX) of an area of ps 
1 degree diameter centered on IRC+10216, which shows 
the object's "astrosphere" (Fig. 0]). It exhibits a textbook 
picture of phenomena expected from the interaction of a 
circumstellar outflow with the ambient interstellar medium: 
an astrosheath, astropa use, a term ination shock, and an 
astrotail (see Fig. 2 of lUetal [2008. for an illustration of 
these terms). An arc-like segment of the termination shock 
is also detected. The arrow in (Fig. shows the direction 
of IRC+10216's proper motion (corrected for the Solar 
motion) and the distance traveled by the star in 5 000 yr. 

Remarkably, the termination shock also manifests itself 
as a dust shell that has been imaged with Herschel at 160 
and 250 /im with the Photodetector Array Camera and 
Spectrometer (PACS) and the SPectral and P hotometric 
Imaging RE c eiver (SPIRE) (jLadial et al.ll2010D . Moreover, 
iDecin et all (j201llh using PACS at 70, 100, and 160 fmi, 
find multiple concentric dust shells with the outermost one 
reaching almost to the termination shock. 

For modeling such interactions between stellar winds 
and the interstellar medium (ISM), one critical parameter 
is the velocity of the star relative to the ISM, u*, which we 
can estimate from the proper motion we have derived in 
Sect. 13.2.21 see Fig. 0] In the FUV image the astrosheath 
forms an incomplete edge brightened shell (the termination 
shock) around the stellar position. However, this shell is 
broken up toward the WSW-SW, in which direction (most 
of) the structures making up the astrotail extend beyond 
its boundaries. We note that this is the direction expected 
for the astrotail for the proper motion vector we determine, 
which has a position a ngle of 70° (E of N) . 

Theoretical studies (Ivan Buren fc McC rav 1988; Wilkin 
|1996|) have yielded a formula for the scale, Zi, at which the 
ram pressure of the freely flowing wind from a moving star 
equals that of the interstellar medium. This is measured as 
the distance between the termination shock and the star 
along the astropause's symmetry axis and is a function 
of mass loss rate, M*, terminal wind velocity, v w , mean 
molecular mass per atom, /2h, ISM number density, nisM, 
and v*. Solving for the stellar velocity, we have 

u* cx (M v w /p, ni SM ) 1/2 /li 

Using the h determined bv lSahai &; Chronopoulosl ((2010) , 
adopting M* of 2 • 10~ 5 M yr -1 (see [T|), and v w 
= 14 km s -1 ), and assuming a low density atomic 
ISM (/2h and n = 1.33 and 1 cm -3 ), respectively, 
ISahai fc Chronopoulosl (|2010| ) derive a stellar velocity of 




Fig. 4. IR C+10216 as seen by GALEX i n 2008 in the FUV 
band (see ISahai fc Chronopoulosl l2010l ). The image has 
been smoothed from the original resolution with a Gaussian 
kernel with a FWHM of 5 pixels. The star symbol indicates 
the position of IRC+10216 and the arrow represents its 
proper motion (corrected for the solar motion) for 5,000 

yr- 



91 km s _1 . lLadial et al.l (|2010f) obtain a slightly higher 
value of 107+9 km s _1 (for «ism = 1 cm -3 ). These values 
assume that the astropause's symmetry axis lies in the sky 
plane. 

To compare the above velocities with our proper motion, 
we must transform our heliocentric value into the galactic 
reference (LSR) frame and correct for the solar motion . 
Following the instructions of lJohnson fc Soderbloml (|1987l ). 
we transform our heliocentric proper motions and radial 
velocity, vhcEL 01 — 19 km s -1 to a velocity vector in 
the galactic frame of reference, (U, V, IU)h c i, which has 
to be corrected for the solar motion, to, finally, obtain 
(U 7 V, W)lsr in the LSR frame. Here, the components U, V, 
and W are positive in the directions of the Galactic center, 
Galactic rotation and the North Galactic pole. If we do this 
transformation, we obtain (U, V, W)lsr — (+34, +23, +10) 
km s _1 for IRC+10216, which has an absolute value of 
42 km s -1 with an estimated uncertainty of less than 10 
km s . 

We now have determined the three dimensional motion 
vector for IRC+10126 in the LSR frame. What we do not 
know is the motion of the ISM in that frame. For the radial 
component of this motion, we can get some insight into this 

We calculate this heliocentric velocity for IRC+10216 
assuming an LSR radial velocity, ?jlsr, of —26 km s~ , using 
the solar motion ve ctor (+10, +5, +7) km s~ , derived from 
Hipparchos data by iDehnen fc Binnevl (|1998T ). Note that our 
value of vlsr differs from the ofte n quoted —22 km s -1 (as, e.g., 
given in the catalog compiled bv lLoup et al.l [19931 ). The latter 
value is derived from low-J CO line profiles that over estimate 
the LSR velocity because of self absorption in their blue wings. 
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question by examining our knowledge of the ISM around 
the star. As discussed in Appendix|Bj the radial component 
of IRC+10216's motion is between 13 and 26 kms" 1 lower 
than that to its ambient ISM. 

Having no other constraints, we assume that this ISM 
is partaking in Galactic rotation with no peculiar velocity 
contribution, i.e., (U, V, WOlsr = (°A°) km s" 1 . Then, our 
best estimate value of the motion of IRC+10216 relative 
to its ambient ISM, 42 km s" 1 , is a factor of sa 2-3 
smaller than the values quote d bv lLadial et al.l (|2010D and 
iSahai fe Chronopoulosl (l2(H0h . 

We note that several assumptions made by the above 
authors could have led to significant uncertainties in their 
estimate of u*. For one, mass loss rates of the order 
of IRC+10216's, 2 • 10" 5 M Q yr" 1 , (or even higher) 
are char acteristic with the latest stages of a s tar's AGB 
lifetime (|Volk et al.ll2000HPerinotto et al.ll2004D . It is thus 
very unlikely, that IRC+10216's mass loss rate had its 
present very high value over the several times 10 4 yr time 
scale it took the astrosphere to be established. In fact, 
co ncentric shell s inCSEs as observed around IRC+10216 
by iDecin et al.l (|201lt) , have been interpreted as evidence 
for variable mas s loss and are p ossibly the result of multip le 
thermal pulses (|01ofsson et al.lll990t IZiilstra et al.lll992f) . 

Second, the a mbient 

ISM's density, for which ISahai fc Chronopoulosl (|2010l ) 
assume 1 cm -3 , is poorly constrained. In the case of Mira 
(o Ceti), which a lso excites FUV em ission while moving 
through the ISM (iMartin et al. 20071) . modeling suggests a 
50 times smaller value ([Wareing "etall[2007l lUetaT [20081. 
Excitation of molecular hydrogen (H2) by hot electrons 
is considered the best candidate mechanism for the FUV 
emiss i on from IRC+10216's astro sphere ()Martin et al.l 
l2007t ISahai fc Chronopoulosl l2010D . The presence of 
molecular hydrogen would imply that IRC+10216 is 
moving within a diffuse molecular cloud. Diffuse molecular 
clouds, ho wever, have hydrogen number densities of 100- 
500 cm -1 (|Snow fc Mc Call 200 61), significantly larger than 
the 1 cm -1 assumed by ISahai fc Chronopoulosl (|2010f ) (and 
1.33 < /2h < 2.33). As a caveat, we point out that Mira's 
FUV emission distribution looks distinctly different from 
IRC+10216's. It has a very long, much more turbulent 
tail with an elongated cometary shape and, while it has 
a pronounced bow shock, it lacks the thin ring making 
up more than a semi circle that marks the boundary of 
IRC+10216's astropause. For the latter, hopefully, further 
modeling will result in a more meaningful estimate of 
the ISM density and other quantities. At present, the 
above facts, plus the unknown inclination angle, prohibit 
a meaningful estimate of v* from the UV observations. 
Vice versa, we note that even if we assumed that v* were 
equal to our estimated speed (42 km s _1 ) the uncertainties 
addressed above would make any conclusions based on the 
UV emission of IRC+10216 and its envelope, e.g., on the 
duration of the AGB phase, extremely uncertain. 

5. Conclusions 

Using the VLA we have imaged the radio emission from 
the archetypical high mass loss carbon-rich AGB star 
IRC+10216. We obtain a precise measure of the size of 
its radio emission (83 mas) and its brightness temperature 
(1660 K). Since the emission is consistent with arising from 
a black body, these quantities, together with the assumed 



distance of 130 pc, determine the average stellar luminosity 
to be 8640 L®. This value is in excellent agreement with 
what is predicted from the period-luminosity relation for 
carbon Miras. 

Our position determination, with ss 10 milli arcsecond 
accuracy, together with older data, allow a determination of 
IRC+10216' s proper motion. After correcting for the solar 
motion, we find that the star moves with a velocity of 39 
km s _1 in the plane of the sky in east-northeastern direction 
(position angle 70°) and has a velocity of 42 km s~ x 
with respect to the LSR. The direction of this motion is 
roughly consistent with FUV and FIR images that show 
an extended astrosphere around IRC+10216, which has a 
broken shell and tail-like features in the opposite direction. 
Our calculation of the star's three dimensional velocity and 
an analysis of the kinematics of its surrounding ISM suggest 
a lower relative velocity than derived in recent studies. This 
suggests a lower (time-averaged) mass loss rate and/or a 
higher ISM density than previously assumed. 
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Appendix A: Phase determination for the epoch of 
the VLA observations 

In order to derive the phase of IRC+10216's light curve 
at the epoch of our VLA observations (2006 February 
23), w e used infrared photometry from IShenavrin et al.l 
(|201lD . This compilation consists of 310 measurements 
in the JHKLM bands obtained between 1999 December 
10 and 2008 November 11 (see Fig. IA~T]) . Each of the 
light curves was interpolated by a cosine fit using an 
implementation of the nonlinear least-squares Marquardt- 
Levenberg algorithm in gnuploiQ- The period (P), the 
epoch of the chosen maximum (JD max ), the average 
magnitude, and the amplitude of variations were treated 
as free parameters. For the HKLM bands, this procedure 
resulted in periods between 628 and 632 days with typical 
uncertainty of 1 day, and the moment of maximum between 
2454550 and 2454560 with typical uncertainty of 4 days. 
The J light curve has slightly poorer quality data and 
the fit gave results slightly different than for the other 
bands (P=638.5 and JD max =2454585.0) and was discarded 
from the further analysis. The combined values for the 
HKLM light curves, in terms of weighed-mean and rms, 
are P=630.0±2.9 days and JD may =24 54554.0±7 4 This 
is consistent with the range of values iLe Bertrel (jl992) 
determined for P from monitoring in different IR bands, 
635-670 d (with an average of 649 d) . His data were taken 
in the mid 1980s over a time range of duration less than 
half than that of our dataset. We next used these values to 
compute the phase of our VLA observations, which gave 
</>=0.79±0.06 (where </>=0 corresponds to the maximum 
brightness). This result agrees very well with the value 
of (^=0. 72 computed from the stel lar luminosity variations 
given in lMen'shchikov et al.l (|2001l) (their Eq. 1), which was 
based on IR variability of the source between 1965 and 
1998. 



6 http://www.gnuplot.info 
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Fig. A.l. Bottom to top: J, H,K,L, and M infrared light 
curves of IRC+10216 (from IShenavrin etail (|2011h ) and 
the our cosine fits to the variability curves. The vertical 
line represent the epoch of the VLA observations (2006 
February 23). 
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Fig. B.l. The profile of the 21 cm line of H I averaged over 
a large region around IRC+10216 chosen to represent the 
typical line of sight emission of neutral hydrogen within a 
few deg rees from the star . Data extracted from the GALFA 
survey (|Peek et alJl201lh . 

Appendix B: On the radial component of the 
motion of IRC+10216 relative to its local ISM 

Insight into the interstellar vicinity of IRC+10216 can 
be gained from HI observations toward the source and 
around it. For this purp ose we extracted a data cube 
from the GALFA survey (|Peek et al.ll201ll) . which covers 
the full astropause of IRC+10216 and a large region 
around it at a resolution of 4'. The HI emission, although 
inhomogeneous on angular scales comparable to the size 
of the astropause, is well represented by the profile shown 
in Fig. IB.ll It is dominated by two components centered 
around -5 kms -1 and +3 kms -1 and with FWHM line 
widths of 11 and 1.5 kms -1 , respectively. The broader 
component has very extended wings seen from -60 to 60 
km s -1 . At some positions around the astropause, a weaker 
component can be seen centered at about -25 kms -1 . There 
are no morphological arguments for any connection of the 
observed emission with the astrosphere of IRC+10216. The 
H I data itself do not allow to distinguish which components 
can be physically collocated with the star, if any. 

Much has been learned about the ISM toward 
IRC+10216 by studies of a bsorption lines seen in optical 
spectra of nearby field stars. [Kendall et all ((2002) identified 
two stars at angular distances of 37" and 137" from 



the position of IRC+10216, which are located behind its 
circumstellar envelope (at spectro/photometric distances 
of 0.5 kpc and 1.4 kpc). The resonance lines of Nal and 
KI show clear signatures of the CSE material but do 
not show any extra component that could be identified 
with the ISM. It cannot be ruled out, however, that some 
atomic absorption is hiding within the strong profiles of the 
circumstellar absorption. The same authors find the diffuse 
interstellar band (DIB) at 6284 A in the spectra of the 
background stars at WHci = —2.9 kms -1 (corresponding 
to wlsr = —9.6 kms -1 ) which they interpret as an ISM 
feature. They find the same velocity component in the 
profile of the Nal lines of a star which is seen outside the 
astrosphere of IRC+10216 (2?5 away from the the stellar 
position) and is located at a distance of 55 pc. This suggests 
that the DIB originates in a nearby, foreground cloud. 
In another approac h to th e same observational material, 
iMauron fc Hugginsl ((2010) found absorption features of 
Call centered at the expected velocity of the circumstellar 
material but broader than the K I lines. One interpretation 
is that interstellar material is partially responsible for those 
features, but a fully circumstellar origin cannot be excluded 
and is actually suggested by the symmetric broadening of 
the profiles. 

These optical studies indicate that the HI emission at 
positive LSR velocities most likely arises in far, background 
clouds, while the component at -10 kms -1 is probably 
related to a foreground cloud. Given the monotonic 
character of Galactic rotation, it is very unlikely that gas 
at even more negative velocities is situated between these 
two mentioned components, so it is reasonable to conclude 
that also clouds at velocities lower than -10 kms -1 are 
foreground objects. The monotonicity of Galactic rotation 
and our results for IRC+10216's motion discussed above 
also imply that the star is moving though a medium 
which has an LSR velocity somewhere between —10 kms -1 
and 3 kms -1 . Given that IRC+10216's LSR velocity is 
—26 kms -1 , the radial component of its motion is between 
16 and 29 kms -1 smaller than that to its ambient ISM. 
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